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ABSTRACT 

We present new H I 21 cm absorption data and ultraviolet spectroscopy from 
HST/STIS of the QSO/galaxy pair 3C 232/NGC 3067. The QSO sightline lies near 
the minor axis and lf8 (llh^ kpc) ab ove the plane of NGC 3067, a nearby luminous 
(cz = 1465 ± 5 km s , L = 0.5 L*) starburst galaxy with a moderate star formation 



rate of 1.4 M yr _1 . The UV spectra show that the Si IV and C IV doublets have the 
same three velocity components at cz = 1369, 1417, and 1530 km s _1 found in Ca II 
0^ ' H & K, Na I D, Mg I, Mg II, and Fe II, implying that the low and high ionization 

gas are both found in three distinct absorbing clouds (only the strongest component 
O ■ at 1420 km s _1 is detected in H I 21 cm). The new Lya observation allow the first 

measurements of the spin and kinetic temperatures of halo gas: T s = 435 ± 140 K and 
Tjj/Tg ~ 1. However, while a standard photoionization model can explain the low ions, 
the C IV and Si IV are explained more easily as collisionally-ionized boundary layers of 
the photoionized clouds. Due to their small inferred space velocity offsets (Av = —260, 
— 130, and +170 km s^ 1 ) relative to the nucleus of NGC 3067 and the spatial coinci- 
dence of low and high ionization gas, we propose that these absorbers are analogous 
^ ■ to Galactic high velocity clouds (HVCs). A comparison of the NGC 3067 clouds and 

Galactic HVCs finds similar H I column densities, kinematics, metallicities, spin tem- 
peratures, and inferred sizes. We find no compelling evidence that any halo gas along 
this sightline is escaping the gravitational potential of NGC 3067, despite its modest 
starburst. 

Subject headings: galaxies: halos — galaxies: individual (NGC 3067) — galaxies: star- 
bursts — quasars: absorption lines — quasars: individual (3C 232) 
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1. Introduction 

Starburst-driven superwinds, powered by the cumulative effect of coeval supernova explosions, 
have long been suggested as a mechanism by which metals and energy can be transported from 
galaxies to enrich their surroundings (e.g., De Young 1978; Heckman 1998). However, observations 
have yet to prove that starburst galaxies commonly generate winds that are capable of escaping the 
galaxy's gravitational potential. Most studies measure the wind velocity near the disk of a starburst 
galaxy. Using such velocities to try to determine whether the wind is gravitationally bound to the 
galaxy is not straightforward, however, since the expanding superbubble will inevitably sweep up 
gas as it travels through the galactic halo, causing the wind to decelerate. 

The spatial extent of wind gas can be measured by Ha or X-ray imaging (Watson, Stanger, 
& Griffiths 1984; Martin 1999; Martin, Kobulnicky, & Heckman 2002) and follow-up spectroscopy 
can determine the wind velocity. The X-ray images themselves allow one to measure the wind 
temperature, which is found to be nearly constant as a function of a galaxy's maximum H I rotation 
speed (Martin 1999), suggesting that superwinds are more likely to escape from dwarf galaxies than 
more massive ones. These emission line studies are limited to studying the densest regions of wind 
gas, however, and cannot measure the full extent of the gas or the wind velocity in the diffuse halo. 

Given a suitably bright background source, typically the stellar continuum of the starburst 
region itself, high velocity wind gas can also be detected in absorption (e.g., Heckman et al. 2000, 
2001). Such studies detect typical outflow velocities of 400-1000 km s _1 , although the outflow 
velocities in dwarf galaxies can be considerably lower (Heckman et al. 2000; Martin 2003). While 
absorption line studies are more sensitive to diffuse gas than emission line studies, they suffer from 
an ambiguity in the distance between the background source and the absorbing gas. These studies 
can only definitively state that the absorber lies in front of the background source and yield no 
information about how far in front of the source the absorber lies. Since absorption line studies of 
galactic winds typically use the starburst itself as a background source, they cannot differentiate 
between a high velocity outflow in the galactic halo or within the starburst region. 

Currently, the only method that directly measures the properties of the superwind gas well 
away from its parent galaxy relies on searching for absorption from the wind gas against a bright 
background QSO (e.g., Stocke et al. 2004). Unfortunately, coincidences of a bright (my < 17) QSO 
and a low-z galaxy separated by < 150hjQ kpc in projection are exceedingly rare; only ~ 60 are 
known to exist over the entire sky, and most of the galaxies in these pairs are not starbursts. Thus, 
a systematic study of galaxy winds with this technique is not yet possible, although much can be 
learned by observing those QSO/galaxy pairs that do exist. 

Rest-frame ultraviolet absorption line spectroscopy of bright QSOs is a very sensitive tracer of 
diffuse gas in interstellar and intergalactic space. For example, the weakest Lya absorption lines 
detected by the Hubble Space Telescope (HST) probe gas at H I column densities, Nm ~ 10 13 cm~ 2 , 
comparable to the weakest Lya lines detected in Keck and VLT ground-based spectra at high 
redshift. Associated metal- line absorptions are commonly found in individual H I absorbers down 



-3- 



to Nm ~ 10 14 cm -2 and pixel addition techniques cumulatively detect Si III, C IV, and O VI 
absorptions down to Nm ~ 10 13 cm~ 2 (Schaye et al. 2003; Aguirre et al. 2004; Simcoe, Sargent, &; 
Rauch 2004). Given that the number density of Lya absorbers greatly increases both with increasing 
redshift and with decreasing column density, a "unified" model for these absorbers places them all 
in the extended halos of bright galaxies (e.g., Steidel 1995; Lanzetta et al. 1995). However, the 
radii of such halos become implausibly large at z ~ 2 for the weak (Nm ~ 10 14 cm -2 ) Lya and 
metal line systems: > 150h^Q kpc for L* galaxies. A more recent interpretation (Stocke, Shull, & 
Penton 2005) suggests that, while the higher column density absorbers could be the bound halos of 
bright galaxies, the weaker absorbers arise in unbound galaxy winds produced by the much more 
numerous low luminosity galaxies (Mb > —16). 

QSO absorption line systems with column densities in the range JV H i = 10 17 - 3 to 10 20 - 3 cm , 
known as Lyman limit systems (LLSs), are associated with bright (L > 0. 1-0.3 L*) galaxy halos as 
shown by Steidel (1995, 1998). Steidel found that 55 of the 58 LLS absorbers with Mg II equivalent 
widths > 0.3 A in his sample at 0.2 < z < 1 were associated with bright galaxies at the same redshift 
within a projected distance of 50hjQ kpc. Follow-up spectra of the underlying stellar kinematics 
in these galaxies (Steidel et al. 2002) and HST imaging and UV spectra of the C IV absorption 
in these systems (Churchill et al. 2000) further indicate that LLSs are the bound gaseous halos of 
luminous galaxies and that in many cases these halos share the kinematics of the underlying stellar 
disk. 

3C 232/NGC 3067 is a well-studied QSO/galaxy pair in which the QSO sightline probes the 
gaseous galaxy halo of a low-z spiral galaxy (Burbidge et al. 1971). The 3C 232 sightline lies near 
the minor axis of NGC 3067 and is located 1'.8 (ll/ifg 1 kpc projected) above the galactic plane. 
NGC 3067 is a luminous (M B ~ -19, L « 0.5 L*; de Vaucouleurs et al. 1991; Marzke, Huchra, 
& Geller 1994) nearly edge-on spiral (i = 68^3 degrees; Rubin, Thonnard, & Ford 1982) with a 
radial velocity of cz = 1465 ± 5 km s~ x (Carilli & van Gorkom 1992, hereafter CvG) and a far 
infrared luminosity of Lfir = 2.4 x 10 43 ergs s~ 4 (Shapley, Fabbiano, & Eskridge 2001), which 
corresponds to a star formation rate (SFR) of 1.4 M yr" 1 (Misiriotis et al. 2004). The SFR of 
NGC 3067 is somewhat less than that of other nearby luminous (L ps L*) starbursts like NGC 253 
(SFR = 3.7 M yr" 1 ), NGC 3079 (SFR = 6.4 M yr" 1 ), and NGC 4945 (SFR = 3.1 M yr" 1 ; 
de Vaucouleurs et al. 1991; Marzke et al. 1994; Shapley et al. 2001; Misiriotis et al. 2004; Grimm, 
Gilfanov, & Sunyaev 2003). Thus, one would expect NGC 3067 to produce a wind similar to but 
perhaps weaker than those of other starburst galaxies in the local Universe. 

H I 21 cm absorption against 3C 232 was first detected near the redshift of NGC 3067 at 
cz = 1418 ± 2 km s _1 by Haschick & Burke (1975), who conjectured that the absorption was 
caused by H I in the outermost regions of NGC 3067. CvG observed 3C 232/NGC 3067 in H I 
21 cm emission at 50" resolution using the VLA and found H I distributed slightly asymmetrically 
(the 21 cm emission extends further to the east of the kinematic center than to the west, and the 
kinematic center is offset from the optical center) at velocities ranging from 1310 to 1635 km s _1 . 
CvG also found a region of extended emission north and south of NGC 3067 with velocities ranging 
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from 1350 to 1550 km s -1 , from which the 21 cm absorption seen toward 3C 232 arises. CvG 
speculated that the extended emission is caused either by extended spiral arms, which may be the 
result of a past encounter with an as-yet undetected companion, or by viewing gas in a polar ring 
that is seen nearly face-on. In the former case we would expect the H I seen above and below the 
plane of NGC 3067 to be rotating in the same direction as the underlying disk; in the latter case 
the space velocity field is less well-constrained. 

Metal lines have also been detected along the 3C 232 sightline at the same velocity as the 
H I 21 cm absorption found by Haschick & Burke. Boksenberg & Sargent (1978) found Ca II H 
and K at cz = 1406 ± 11 km s _1 with the Hale 5m. Stocke et al. (1991, hereafter S91) used the 
Multiple Mirror Telescope to find three separate velocity components in Ca II H & K and Na I D 
at cz = 1369 ± 2, 1417 ± 2, and 1530 ± 10 km s" 1 and Tumlinson et al. (1999, hereafter T99) used 
the Goddard High Resolution Spectrograph (GHRS) aboard HST to discover that Mg I, Mg II, 
and Fe II absorption lines show evidence for the same three velocity components that S91 found in 
Na I and Ca II, although the Mg II lines are highly saturated and could cover this entire velocity 
range. None of these spectra preclude weaker absorption at velocities intermediate to the three 
strong systems. 

The NGC 3067 halo cloud probed by the 3C 232 sightline has been used to set both an upper 
limit based upon Ha imaging (I„ < 4 x 10~ 23 ergs s _1 cm~ 2 Hz -1 sr _1 ; S91) and a lower limit 
based upon Mg Il/Mg I absorption line ratios il u > 1 x 10~ 23 ergs s _1 cm" " 2 Hz" 1 sr^ 1 ; T99) on 
the extragalactic ionizing flux (Ij,). The Ha images also yielded a limit on the escape fraction of 
ionizing radiation from the disk of NGC 3067 to the halo locations probed by the 3C 232 sightline 
of fesc < 2% (T99). This value is comparable to the most recently derived f esc = 1-2% for the 
Galaxy from Ha images of high velocity clouds (HVCs; Putman et al. 2003). 

This paper presents H I 21 cm and ultraviolet spectra of 3C 232 to determine if the multiple 
velocity components found by S91 and T99 along this sightline are present in H I Lya, Si IV, and 
C IV. In §2 we describe our radio and UV observations and analysis. In §3 we hypothesize that the 
multiple velocity components found by S91 correspond to HVCs around NGC 3067 and compare 
the properties of these clouds to those of Galactic HVCs. In §4 we summarize our findings and 
discuss the evidence for unbound starburst winds in light of our results. 

2. Observations and Data Reduction 

Our full dataset includes previously-published optical and UV spectra of 3C 232 (S91; T99) 
and published Very Large Array (VLA) H I emission maps (CvG), as well as new H I 21 cm 
observations with the Very Long Baseline Array (VLBA) and the Arecibo radio telescope and 
new HST low and medium resolution far-UV spectra obtained with the Space Telescope Imaging 
Spectrograph (STIS). The H I 21 cm and HST/STIS observations are described below, while the 
previously published data will be described as needed in §3. 
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2.1. VLB A Observations 

High angular resolution observations of 3C 232 were carried out at 1413.7 MHz on 2002 January 
31 using the VLBA, which is operated by NRAO 5 . The bandwidth of the observations was 8 MHz, 
in each of the right- and left-hand circular polarizations, sampled at two bits and centered at the 
frequency of the H I 21 cm line, at a heliocentric redshift of z = 0.004736, or cz = 1420 km s -1 . 
The data were correlated at the VLBA correlator in Socorro, New Mexico, with 1024-point spec- 
tral resolution per baseband channel and 4 s correlator integration time. The resulting channel 
separation was 7.8 kHz (1.7 km s -1 ). The total observing time was 11.7 hr. Along with 3C 232, 
the radio source J0927+3902 was observed at the same frequency as a fringe finder and a bandpass 
calibrator. Data reduction and analysis were performed using the Astronomical Image Processing 
System (AIPS), and the Astronomical Information Processing System (AIPS++), of NRAO. 

Figure 1 shows our high resolution L-band continuum image of 3C 232 with 10.7 x 6.7 mas 
resolution, along with H I absorption spectra averaged over two regions. The image was obtained 
by using the CLEAN algorithm as implemented in AIPS, with an intermediate grid weighting 
(Robust=+2 in IMAGR). The continuum source consists of a compact core and a jet structure 
extending to the south, with a weaker extension to the east (but see Figure 2). A small amount 
of additional, very low surface brightness flux exists on the 500 mas scale (Condon et al. 1998). 
Table 1 lists the results of fitting Gaussian models to the observed spatial profile for the source 
components listed in column 1. We designate the two regions of 3C 232 toward which the spectra in 
Figure 1 have been extracted as core and jet. No H I absorption larger than 3a = 4.7 mJy beam -1 
is detected toward other continuum regions in 3C 232. The velocity full width at half maximum 
(FWHM) of the H I absorption is 4.7 ±0.1 km s -1 centered at 1420.0 ±0.1 km s -1 . No H I 
absorption is detected at 1369 and 1530 km s -1 , in contrast to where metal lines appear in optical 
and UV data (S91; T99). 

Figure 2 is a wide field continuum image with the same angular resolution as Figure 1, showing 
an extended structure at an angular distance of ~ 180 mas to the east of the main continuum 
structure of 3C 232. A gaussian fit to this structure results in a total flux density of 46.9 ±0.4 mJy, 
and a size at full width half maximum of 60 x 36 mas on the plane of the sky. No H I absorption 
larger than 3a = 4.7 mJy beam -1 is detected toward this diffuse continuum emission region. 

The peak optical depth values in the core and the jet regions are 3.5 ± 0.8% and 2.6 ± 0.8%, 
respectively, which yields column densities (Ahi) of H I absorption of Nm/T s = 1.7(±0.4) x 
10 17 cm -2 K -1 and 1.1(±0.3) x 10 17 cm -2 K -1 for the core and jet regions, respectively. An 
absorber with this column density would cause a flux decrement of ~ 0.1 mJy beam -1 against the 
eastern extension of the jet (see Figure 1) and ~ 1 mJy beam -1 against the diffuse emission region 
in Figure 2. Since the 3a detection limit in these regions is 4.7 mJy beam -1 , we cannot exclude the 



5 The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation operated 
by Associated Universities, Inc. 
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Fig. 1. — HI absorption spectra obtained at two locations against the background continuum 
source 3C 232 at 1413.7 MHz. No absorption > 3a is detected toward other continuum emission 
regions in 3C 232. The restoring beam size is 10.7 x 6.7 mas in position angle —0.2°. The rms 
noise level in the H I spectra is 1.6 mJy beam -1 . The velocity resolution is 1.7 km s _1 . The 
contour levels of the continuum image are —5, 5, 10, 20, 40, 2560 times the rms noise level, 
which is 123 fiJy beam -1 . The reference position (0,0) is a(J2000) = 09 h 58 m 20?9499, <5(J2000) = 
32° 24' 02'/177. At the redshift of NGC 3067, 1 mas = 0.1 pc for H = 70 km s -1 Mpc" 1 . 

possibility that the H I absorption is constant over a region extending at least 20 mas (2/iyg 1 pc) 
and possibly as large as 0.2" (20h^ pc) on the sky. 

CvG measured the column density of H I 21 cm emission toward 3C 232 to be 8(±4) x 10 19 cm -2 
using the VLA (50" beam). Comparing this value with the Nm/T s derived from our VLBA H I 
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Fig. 2. — A large-field continuum image of 3C 232 at 1413.7 MHz. The restoring beam and contour 
levels are the same as in Figure 1 and the peak flux is 560 mJy beam -1 . Note the diffuse continuum 
emission ~ 180 mas to the east of 3C 232. 

absorption observations, we derive a hydrogen spin temperature of T s = 470 ± 260 K. Of course, 
this value assumes that the H I 21 cm emission and absorption have comparable physical extent. 

2.2. Arecibo Observations 

Single dish observations of 3C 232/NGC 3067 were carried out with the L-band Wide receiver 
of the 305 m Arecibo radio telescope 6 on 2003 November 23-24, for a total of 2.5 hr. The double 
position switching observations utilized two of the interim correlator boards to observe at the 
frequency of the H I 21 cm line in NGC 3067. The bandwidth of each board was 6.25 MHz, with 
one linear polarization each and 2048 spectral channels. The source 4C+32.34 was observed as a 
bandpass calibrator. 

Figure 3 is a total intensity Hanning-smoothed H I spectrum of 3C 232/NGC 3067. The 
plot clearly shows the 300 km s _1 wide H I 21 cm emission in NGC 3067, and the narrow 



The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell 
University under a cooperative agreement with the National Science Foundation. 
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Fig. 3. — Hanning-smoothed H I 21 cm spectrum of the galaxy NGC 3067 taken at Arecibo. The 
velocity resolution is 1.3 km s _1 , and the rms noise level is 1.2 mJy The arrows indicate the 
positions of the other two velocity components (cz = 1369, 1530 km s^ 1 ) found in various metal 
species (S91; T99). Note that all three radial velocities fall well within the H I emission profile of 
NGC 3067, and that the H I 21 cm absorption detection is close to the galactocentric radial velocity 
of 1465 km s" 1 (CvG). 



H I 21 cm absorption line that arises on the line of sight to the background quasar 3C 232 at 
1420.2 ±0.1 km s _1 . The velocity resolution of this spectrum is 1.3 km s _1 , and the rms noise level 
is 1.2 mJy. The velocity FWHM of the 21 cm absorption is 4.2 ± 0.1 km s _1 and its peak optical 
depth is 2.81 ± 0.08%, yielding a column density of N m /T s = 2.3(±0.1) x 10 17 cm -2 K _1 . While 
the single dish optical depth is marginally larger than the VLB A measurements, these observations 
suggest no strong variation in covering factor or T s over a size scale of 0'/2, or 20/i^q 1 pc. Again 
using the H I 21 cm emission column density found by CvG, this value corresponds to a hydrogen 
spin temperature of T s = 350 ± 175 K, which agrees with the spin temperature calculated from the 
VLBA column density to within errors. 

No H I absorption is detected at the two other metal- line velocities to r < 0.24% (3cr). Assum- 
ing a spin temperature of T s = 450 K and the same velocity width as the 1420 km s _1 absorption 
line, this optical depth limit corresponds to a 3a column density limit of Nm < 9 x 10 18 cm~ 2 . 
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Notice that the three metal line velocities are well within the extent of the disk rotation curve. 

2.3. HST/STIS Observations 

3C 232 was observed by HST/STIS 7 with the G140L grating as part of GO Snapshot Program 
9506 (J. Stocke, PI), and also with the G140M grating at three wavelength settings as part of 
GO Program 8596 (M. Pettini, PI). A journal of these observations appears in Table 2; total 
exposure times were 1.2 ksec with G140L, 11.2 ksec with G140M at central wavelengths of 1222 
and 1400 A, and 22.3 ksec with G140M at a central wavelength of 1567 A. Multiple exposures 
with the same central wavelength were shifted to a heliocentric reference frame and coadded with 
signal-to-noise weighting. The continuum in each wavelength region was normalized using Legendre 
polynomials and the resulting spectra were examined for evidence of absorption from NGC 3067. 
These normalized spectra are shown in Figure 4 and absorption lines associated with NGC 3067 
are labelled at their observed wavelengths. 

The G140L spectrum yields data on Lya, O I, C II, C IV, Si II, Si III, Si IV, S II, S III, 
and Al II. All of the prominent Si and C lines in the wavelength range covered by this spectrum 
are detected over a wide range of ionization states. The resolution of this spectrum is too coarse 
(250 km s" 1 ; see Table 2) to separate the blended S Ill/Si II and S Il/Si II lines; however, since the 
oscillator strengths of the Si II lines in these blends are 1-2 orders of magnitude stronger than those 
of the S lines and we detect no S II A1251 or S II A1254 absorption in our spectrum, these blended 
lines primarily reflect the properties of the Si II ions. Likewise, these values are the combined 
equivalent widths of all three velocity components detected in other metals in this sightline, as the 
resolution of this spectrum is too coarse to separate the components found by S91 and T99, which 
have a maximum separation of Av = 160 ± 10 km s _1 . Therefore, we analyzed the higher resolution 
(20-25 km s _1 ) G140M spectra to search for these velocity components. Measured radial velocities 
and equivalent widths for the species detected in the G140L spectrum are shown in Table 3. 



7 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS5-26555. 
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Table 1. Gaussian Fits to the Continuum Components in 3C 232 (Fig. 1) 





Relative Position 3, 


Peak 


Total 


Major b 


Minor b 


PA. 


Source 


(mas) 


(mJy beam -1 ) 


(mJy) 


(mas) 


(mas) 


(°) 


core 


, 


478.2 ± 0.3 


587.0 ±0.3 


8.7 


4.7 


175 ±5 


jet 


1.5W, US 


196.2 ±0.5 


388.2 ±0.9 


11.0 


6.1 


179 ±5 


east. ext. 


7.5E, 14S 


3.92 ± 0.02 


4.66 ± 0.03 


8.3 


4.8 


43 ±20 



a The position (0,0) is a(J2000) = 09 h 58 m 20?9499, 5(J2000) = ±32° 24' 02T177. 
b Errors are « 1 mas 



Table 2. Journal of the HST/STIS Observations Toward 3C 232 



Wavelength Coverage Resolution Usable Exposure Time 



Setting (A) (km s x ) Date of Observation (s) 



1 1372-1428 20 2001 Oct 29 11,178 

2 1194-1250 25 2001 Oct 30-31 11,178 

3 1540-1594 20 2002 Jan 02 8,208 

1540-1594 20 2002 Jan 04 14,098 

4 1150-1700 250 2003 Mar 04 1,200 
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Table 3. Best-Fit Parameters to the Lines Detected in the STIS Spectra 



Ion 






log N h 


w A 




(A) 


(km s _1 ) 


(cm -2 ) 


(A) 


S Ill/Si II . 


1190.2, 1190.4 


2230 ± 100 




2.0 ±0.3 


Si TTT 


1 206 5 


1 ci n _i_ i n 


17 a _i_ n i 


1 7 + 3 


Lya 


1215.7 


1430 ± 20 


20.0 ±0.1 


6.5 ±2.4 


S Il/Si II . . 


1259.5, 1260.4 


1550 ± 120 




1.2 ±0.2 


01 


1302.2 


1520 ± 110 




0.4 ±0.1 


Si II 


1304.4 


1350 ± 210 




1.6 ±0.3 


C TT 


1 334 5 


1 ' ) >J \ I _l_ \J\J 




1.4 ± 0.1 


CI • nr 

Si IV 


1393.8, 1402.8 


1 ot r i on 

1375 ± 20 


13.0 ± 0.9 


0.1 ± 0.1 






1435 ± 15 


16.8 ±0.5 


0.4 ±0.2 






1535 ± 10 


15.8 ± 1.5 


0.4 ±0.2 


Si IV Total: 






16.8 ±0.8 


0.9 ±0.3 


Si II 


1526.7 


1580 ± 50 




1.3 ±0.4 


C IV 


1548.2, 1550.8 


1375 ± 20 


16.2 ±2.2 


0.3 ±0.2 






1455 ± 20 


18.2 ±0.5 


1.1 ±0.6 






1560 ± 20 


15.9 ±2.3 


0.4 ±0.2 


C IV Total: 






18.2 ±0.8 


1.8 ±0.7 


Al II 


1670.8 


1430 ± 90 




1.2 ±0.3 



a Tabulated recession velocities are heliocentric. 

b Lines without a reported column density were measured from 
the G140L spectrum. All other lines were measured from the G140M 
spectra. 
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Fig. 4. — The final coadded, normalized STIS spectra. The top panel shows the G140L snapshot 
spectrum, the second panel shows the G140M spectrum with a central wavelength of 1222 A, the 
third panel shows the G140M spectrum centered at 1400 A, and the bottom panel shows the G140M 
spectrum centered at 1567 A. Absorption features from NGC 3067 are labelled in each panel; in 
the G140L spectrum, the Si III and Lya absorption from NGC 3067 are blended with Galactic Lya 
absorption. In the G140M spectra, features marked with a "G" are absorption from our Galaxy 
and those marked with an "i" are absorption intrinsic to 3C 232 at z = 0.533. Features marked 
with a "?" are likely low-z Lya absorbers since the Lya forest region of 3C 232 extends to 1864 A. 
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Voigt profiles were fitted to the Lya, Si III, Si IV, and C IV absorption features in the G140M 
data. While the Lya absorption shows no sign of velocity structure, the small velocity offsets 
associated with the multiple metal-line components would not be detectable due to the wide damped 
Lya profile. No clear velocity substructure is present in the Si III A1207 absorption feature from 
NGC 3067 at A ps 1213 A either, although the trough of the absorption feature spans velocities 
from cz ps 1400-1650 km s" 1 (see Figure 5). The expected location of the N I A1200 absorption 
feature at the redshift of NGC 3067 lies in the middle of the Galactic Si III A1207 absorption, so 
we cannot comment on its strength. 

In order to determine the most robust values of the recession velocity and the column density, 
the observed Lya profile was fit by a Voigt profile with constant 6- value and varying cz and log N 
in order to minimize the x 2 °f t ne fit. The equivalent width of the best-fit profile was calculated 
from the corresponding model parameters, not from a direct non-parametric measurement of the 
absorption. The best-fit recession velocity, column density, and equivalent width for the Lya profile 
are shown in Table 3 and the corresponding Voigt profile is shown in the top-left panel of Figure 5. 
The best-fit recession velocity of the Lya profile agrees with the H I 21 cm velocity to within the 
errors. The best-fit column density (iVm = 10 20 ' cm~ 2 ) makes this absorber a Lyman limit system, 
which are typically associated with the halos of bright galaxies (Steidel 1995, 1998; Stocke et al. 
2005). Together with the H I 21 cm radio data, the Lya column density yields estimates of the spin 
temperature from the single dish and VLBA data of 435 ±140 K and 590 ±380 K, respectively. The 
Si III profile was fit in the same way as the Lya profile and its best-fit recession velocity, column 
density, and equivalent width are shown in Table 3. The corresponding Voigt profile is shown in 
the bottom-left panel of Figure 5. 

In contrast to the H I 21 cm and Lya data, the C IV and Si IV profiles have multiple velocity 
components as seen in several other metal species in the 3C 232 sightline (S91; T99). Multiple 
Voigt profiles were fit to the Si IV and C IV doublets in order to determine the most robust values 
of recession velocity and column density. The initial velocity estimates for these profiles were cz = 
1370, 1420, and 1530 km s _1 (chosen to match the velocity components at cz = 1369 ± 2, 1417 ± 2, 
and 1530 ± 10 km s" 1 found in Na I and Ca II by S91). A range of b- values from 4-20 km s was 
chosen and held constant while the recession velocity and column density were allowed to vary in 
order to minimize the x 2 of the fits to the Si IV and C IV profiles. The best-fit parameters from 
the fits with different 6-values were then compared to estimate the fitting errors. The resulting 
best-fit recession velocities, column densities, and equivalent widths are shown in Table 3. The 
corresponding Voigt profiles for the Si IV and C IV doublets are shown in the middle and right 
panels of Figure 5, respectively. The best-fit velocities for the 1370 and 1530 km s~ x components in 
both the Si IV and C IV profiles agree with the values found by S91 to within the combined errors. 
The best-fit velocity at cz = 1435 ± 15 km s _1 in the Si IV profile is just barely discrepant from 
the S91 value of cz = 1417 ± 2 km s" 1 at the la level, while the best-fit velocity at cz = 1455 ± 20 
in the C IV profile differs from the S91 velocity by 1.8(7. Considering that the Si IV and C IV 
components are significantly blended (although much less so than Mg II; T99), all of the best-fit 
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Fig. 5. — Voigt profiles generated using the recession velocity and column density values from 
Table 3 and overlaid on the STIS G140M spectra. Each panel is labelled with the absorption 
feature that it depicts. The dotted vertical lines show the three velocity components at cz = 1369, 
1417, and 1530 km s _1 found in Ca II H & K and Na I D by S91. Note that the Lya panel shows 
a more extended velocity range than the other panels because the profile is black over nearly all of 
the narrower velocity range. 



velocities agree well with the component velocities determined previously for Ca II and Na I (S91). 



3. Discussion 



Perhaps the most remarkable aspect of this sightline is the consistent velocity structure from 
neutral elements (Na I D, Mg I) to singly-ionized species (Ca II H & K, Mg II, Fe II) to high 
ions (Si IV, C IV). S91 found three velocity components in the profiles of Na I D and Ca II H & 
K and determined that these components lie at heliocentric recession velocities of cz = 1369 ± 2, 
1417 ± 2, and 1530 ± 10 km s _1 . T99 used HST/GHRS to examine the absorption-line profiles 
in the wavelength region A = 2250-3000 A and found absorption lines of Mg I, Mg II, and Fe II 
associated with NGC 3067. The Mg I A2853 profile has velocity components at heliocentric velocities 
of cz = 1370 and 1420 km s" 1 and the Fe II A2344, A2374, and A2600 profiles all show broad 
absorption centered at cz = 1450 km s -1 and extending from cz « 1350-1600 km s _1 , although no 
clear velocity substructure is present (the Fe II A2586 line was also observed, but due to blending 
with Galactic Fe II A2600, its velocity structure is indeterminate). Like the Lya profile described 
above, the Mg II AA2796, 2800 profiles found by T99 are highly saturated. Thus, while absorption 
is present in these profiles from cz = 1300-1600 km s^ 1 , it is impossible to detect, or rule out, 
any velocity substructure in these transitions. Finally, the present study uses HST/STIS spectra 
to examine the Si IV and C IV absorption lines associated with NGC 3067 and finds evidence of 
multiple velocity components in these profiles. We find three distinct velocity components in both 
doublets, and by fitting Voigt profiles to the observed absorption we have found that the velocity 
components for the Si IV and C IV transitions generally match the components found by S91 to 
within the combined errors. We believe, but cannot prove definitively due to line blending, that 
the velocities of Si IV and C IV match the low ion velocities exactly. 

Because all of the far-UV absorption in Table 3 is in the Lya forest region for this QSO 
(z = 0.533), line blending with low-z Lya interlopers (e.g., the line just blueward of Si IV A1394 in 
Figure 4) can contaminate the wavelength region of interest and create false velocity components. 
Lya interlopers are especially problematic for the Si IV and C IV profiles due to the sensitivity of 
the STIS G140M data. Nonetheless, the presence of a consistent velocity structure over numerous 
ionization states implies that the high ionization gas is spatially coincident with the neutral and 
singly-ionized gas; i.e., the gas exists in three "clouds" above the disk of NGC 3067. 

The Na/Ca ratio for these absorbers (iV(Na I)/iV(Ca II) ~ 1; S91) is not consistent with gas 
originating in the Galactic disk (iV(Na I)/JV(Ca II) « 2-200; Crawford, Barlow, & Blades 1989), 
but no studies have measured the Na/Ca ratio out to a Galactocentric radius of 30h^Q kpc, the 
galactocentric radius of the NGC 3067 clouds if they lie in the disk of NGC 3067 rather than its 
halo. However, Morton, York, & Jenkins (1986) found a tentative cutoff in the distance at which 
Ca II could be detected from a galaxy center of ~ D(0)25, the isophotal face-on diameter of a 
galaxy at a surface brightness of 25 mag arcsec" 

~ 2 (D(0) 2 5 = 12/ifo 1 kpc for NGC 3067) and Briggs 
et al. (1980) found that disk gas detected in H I 21 cm at galactocentric radii > 1.5 Holmberg 
radii is uncommon in the local Universe and that the distribution of the disk gas shows a sharp 
cutoff at A^hi < 10 19 cm~ 2 . These results suggest that the NGC 3067 clouds are likely associated 
with the halo (i.e., not the very extended disk) of the galaxy. We propose that each of the velocity 
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components in the 3C 232 sightline is caused by an absorbing cloud of gas above the plane of 
NGC 3067 that is analogous to a Galactic HVC. In §3.1 we compare the kinematics of the absorbing 
clouds around NGC 3067 and Galactic HVCs, in §3.2 we compare their spin temperatures and cloud 
sizes and in §3.3 we discuss the ionization conditions necessary to produce the observed equivalent 
width ratios. 

3.1. Gas Kinematics 

The H I 21 cm emission map and rotation curve of CvG (their Fig. 7) is reproduced in Figure 6, 
where the halo gas to the NE (SW) can be seen to have systematically lower (higher) radial velocities 
than the disk "below" these regions. If the halo gas in these regions is rotating similarly to disk gas, 
then the northern half is closer to us than the southern half (i.e., we are looking at the "underside" 
of the galaxy). Using the systemic velocity for NGC 3067 of 1465 ± 5 km s _1 (CvG), an inclination 
of 68±j de grees (Rubin et al. 1982), and assuming that the clouds are on purely radial trajectories 
along the minor axis of NGC 3067 as one might expect for an outflowing starburst wind, the 1370, 
1420, and 1530 km s" 1 clouds are moving at -260±^, -1301^, and +170±§g km s" 1 with respect 
to the nucleus of NGC 3067, respectively. Thus, given the galaxy's orientation, two of the clouds 
are infalling and one is outflowing. The uncertainties in the above velocities take into account the 
uncertainty in the velocity of each component and the uncertainty in the systemic velocity and 
inclination of NGC 3067, but they do not take into account uncertainties in our physical model. If 
the clouds are moving on radial orbits that are not along the minor axis of NGC 3067 then their 
inflow/outflow velocities could be much larger than this if they are moving nearly along the plane 
of the sky. Of course, these clouds could also share some or all of the rotation of the disk, in which 
case the radial inflow/outflow velocity is diminished. 

The rotation curve in Figure 6 yields a dynamical mass for NGC 3067 of 1.2(±0.2) x 10 u k£ M , 
which leads to an escape velocity of 310±30 km s _1 at a radius of 11 kpc. This velocity implies that 
even if our assumed geometry is wrong and the halo gas is counter-rotating with respect to the disk 
(i.e., the southern side of NGC 3067 is closer to us than the northern side) the H I detected cloud 
at cz = 1420 km s _1 would still be bound to NGC 3067 if it is traveling on a radial orbit along the 
galaxy's minor axis. Especially considering that these clouds are likely to share some rotational 
velocity with the underlying disk, it is likely that all three clouds are bound to the galaxy irrespec- 
tive of the assumed geometry, lending more weight to the supposition that these clouds are not the 
result of an outflowing galactic wind. Despite the modest starburst (1.4 M yr" 1 ) in NGC 3067, 
there is no strong evidence that halo gas is escaping, and since these velocities are shared by both 
the low and high ions, this conclusion is not based solely upon the velocity of the neutral halo gas. 
Note in Figure 6 that the 3C 232 sightline samples a chord close to the minor axis of NGC 3067, 
the expected location of an outflowing starburst wind; e.g., Heckman, Armus, & Miley (1990) find 
that the high velocity gas in several nearby starbursts can be explained by a bi-conic outflow with 
a typical opening angle of ~ 65°. The location of 3C 232 relative to NGC 3067 (~ 10-15° from the 
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Fig. 6. — H I emission map and rotation curve of NGC 3067 reproduced from CvG. The gray 
scale shows the integrated H I column density for NGC 3067 at 50" (5hj^ kpc) resolution, ranging 
from 2.2 x 10 19 atoms cm -2 (white) to 1.3 x 10 21 atoms cm~ 2 (black). The contours show the 
intensity-weighted mean velocity for the H I gas. Contour levels are 1430, 1440, 1560 km s _1 . 
A cross marks the position of 3C 232. 

minor axis; see Figure 6) places this sightline well within this cone angle. 

The usual convention for Galactic HVCs is to call gas with a radial velocity > 90 km s _1 
with respect to the local standard of rest (LSR) an HVC; however, the range of radial velocities 
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allowed by differential galactic rotation varies with Galactic position (Wakker & van Woerden 1997). 
Converting the LSR velocity to the Galactic standard of rest (GSR), an HVC with vlsr = 90 km s" 1 
would have vqsr i n the range -130 to +310 km s" 1 (Wakker & van Woerden 1997). However, 
it should be stressed that this vqsr is a radial velocity only and does not represent the full space 
velocity of the HVC. The observed radial v GS r for the 1370, 1420, and 1530 km s" 1 clouds (-95 
to +65 km s" 1 ) are clearly in this velocity range, and so have kinematics and locations consistent 
with Galactic HVCs. Additionally, many Galactic HVCs (e.g., Complex C; Tripp et al. 2003) 
are observed to be infalling onto our local Galactic disk. Since NGC 3067 also shows evidence 
for infalling clouds and its starburst is of comparable strength to those of other bright (L ~ L*) 
local starburst galaxies, one would expect that their superwinds might not be able to escape their 
gravitational potentials either. 

We have measured a column density of logAni = 20.0 ± 0.1 cm" 2 for the 1417 km s _1 cloud, 
and have calculated an upper limit of Nm < 9 x 10 cm" - 2 for the 1370 and 1530 km s" 1 clouds, 
assuming that they have a similar spin temperature and velocity width as the 1420 km s" 1 cloud. 
Large angle sky surveys such as Bajaja et al. (1985) and Hulsbosch & Wakker (1988) search for 
Galactic HVCs down to a detection limit of Nm ~ 2-3 x 10 18 cm" 2 , so it is questionable whether the 
1370 and 1530 km s" 1 clouds would be detected in such a survey if they were located near the Milky 
Way. The combined Bajaja et al. (1985) and Hulsbosch &: Wakker (1988) surveys find a covering 
factor for Galactic HVCs, as seen from Earth, of ~ 0.4 (Wakker & van Woerden 1991). Thus, if 
NGC 3067 has the same covering factor of HVCs as the Milky Way, there is a ~ 64% probability of 
intercepting at least one cloud with Nm >2x 10 18 cm" 2 along the sightline to 3C 232, assuming 
that the clouds are randomly distributed. Thus the coincidence of one strong HVC and two weak 
ones in this sightline is consistent with the expectations from our Galaxy's HVC population. 

Recently, Thilker et al. (2004) have discovered an extensive population of H I clouds surround- 
ing M 31 with the Green Bank Telescope. They find > 20 discrete clouds within 50 kpc of the 
M 31 disk that they believe are analogs to Galactic HVCs. The M 31 clouds have radial velocities 
comparable to the velocity of the outer disk of M 31 and masses of (0.15-1.3) x 10 6 M (Thilker 
et al. 2004). Additionally, Thilker et al. find a filamentary halo component concentrated near the 
M 31 systemic velocity that extends at least 30 kpc from its nucleus. If the extensive H I cloud 
systems of the Milky Way and M 31 are ubiquitous in the low-z Universe then it is no surprise that 
we have detected similar clouds around NGC 3067. 

3.2. Spin Temperature and Cloud Size 

Our dataset allows us to make several estimates of the spin temperature of the H I in the clouds. 
From our VLBA observations (§2.1) we have calculated a column density for H I absorption toward 
3C 232 of N m /T s = 1.7(±0.4) x 10 17 cm" 2 K" 1 . Combining this result with the column density of 
H I emission (8(±4) x 10 19 cm" 2 ) found by CvG (see Figure 6), we calculate a spin temperature of 
T s = 470 ± 260 K for the strongest H I component seen on the line of sight to 3C 232. By fitting 
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a Voigt profile to the Lya absorption associated with NGC 3067 in our HST/STIS spectrum of 
3C 232 (§2.3), we have found a column density of log A^m = 20.0 ± 0.1 cm -2 for the 1420 km s _1 
cloud, yielding a spin temperature of T s = 590 =L 380 K when combined with the VLBA Nm/T s 
measurement. A third value is obtained by comparing the single-dish Arecibo optical depth with 
the Lya column density: 435 ± 140 K. While this latter value has the smallest statistical errors, the 
radio and UV beam sizes against which the absorption is measured are not as good a match (radio 
continuum size > 30 mas; UV continuum size < 1 mas). In contrast, the VLBA T s value uses only 
the core of 3C 232, which has a size of ~ 10 mas. However, within the rather large measurement 
errors, all of these values agree. 

Unfortunately, very little is known about the spin temperatures of Galactic HVCs, but Akeson 
& Blitz (1999) have measured a 3a lower limit of T s > 180-320 K on the line of sight toward 
one Galactic HVC. The spin temperature measurements of the strongest absorbing cloud near 
NGC 3067 are clearly consistent with this lower limit. On the other hand, H I clouds in the Outer 
Arm of the Galactic disk (radius ~ 15-20 kpc) have spin temperatures that are much lower than 
those calculated for the 1420 km s _1 cloud (typically ~ 50-150 K, c.f., Akeson & Blitz 1999), further 
indicating that this cloud is not associated with the disk of NGC 3067. 

We can also estimate the kinetic temperature (Tk) of the H I using the width of the 21 cm 
absorption profile. This profile has a FWHM of 4.2 ± 0.1 km s _1 in the Arecibo data, which 
corresponds to a b- value of 2.5 ± 0.1 km s -1 and yields Tk = 380 ±30 K. If one performs the 
same calculation starting with the FWHM of the VLBA 21 cm absorption profile then Tk = 
475 ± 35 K. These estimates agree with the spin temperature measurements above to within the 
errors, suggesting that T s and Tk are coupled. To our knowledge, this extragalactic cloud is the 
first for which both T s and Tk have been measured and found to be approximately equal. Particle 
collisions naturally drive the spin temperature of a gas toward the kinetic temperature. Using 
the collisional de-excitation rate coefficients of Allison & Dalgarno (1969), a particle density of 
n Hi ^ 0.1 cm~ 3 is required for collisions to be the dominant source of temperature coupling at 
this kinetic temperature (Kulkarni & Heiles 1988). T99 found a density of nui ~ 0.005 cm -3 for 
the NGC 3067 clouds from photoionization modeling, smaller than the density required for particle 
collisions to couple T s and Tk. 

Resonant scattering of Lya photons can also drive the spin temperature to the kinetic temper- 
ature, since the radiation field exchanges energy with the gas during the scattering process through 
atomic recoil (Field 1958, 1959). This interaction causes T s to approach the color temperature of 
the radiation, which is equal to Tk if enough scatterings have occurred (i.e., the optical depth is 
large enough; Field 1959; Deguchi & Watson 1985). Deguchi <fe Watson (1985) found that Lya 
scattering dominates the excitation of low-density extragalactic clouds if the Lya optical depth 
(t l ) is > 10 5 . Using Equation 7 of Wolfe (1979), r L in the NGC 3067 clouds is 2.9(±0.1) x 10 7 . 
Since the optical depth is so large, the Lya photons cannot penetrate easily into the cloud from an 
external source; thus, the Lya photons must be generated inside the cloud itself from the recom- 
bination of electrons and protons to form H I (Deguchi & Watson 1985; Kulkarni & Heiles 1988). 
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In statistical equilibrium, the recombination rate is equal to the ionization rate of H I by cosmic 
rays and EUV and soft X-ray photons that are capable of penetrating to the cloud's interior. For 
resonant Lya scattering to drive T s to within 10% of Tk at this kinetic temperature, an ionization 
rate of > 5 x 10 -17 ionizations per hydrogen atom per second is required (Deguchi & Watson 1985). 
Watson & Deguchi (1984) found an ionization rate of 2.5 x 10 -18 ionizations per hydrogen atom per 
second for penetrating radiation from the extragalactic 1-10 keV X-ray flux. The true ionization 
rate in the NGC 3067 clouds is likely sufficient to couple T s and T\ (i.e., > 5 x 10~ 17 s^ 1 ) when 
one includes the contributions of cosmic rays, soft X-rays produced in NGC 3067, the extragalactic 
X-ray background, and some Lyman continuum photons escaping from NGC 3067. 

The diversity of our observational data also allows us to estimate the size of the strongest 
absorbing cloud. Since absorption from the 1420 km s _1 cloud is present to comparable optical 
depths in the spectra of both the core and jet regions of 3C 232 (see Figure 1), the absorbing 
cloud must be larger than the separation between the two regions, 20 mas (see Table 1). At the 
redshift of NGC 3067 this angular distance corresponds to a physical distance of 2/i^q 1 pc. Since 
our estimates of T s using either the Lya Nm value or the H I 21 cm emission iV HI value (CvG) 
agree to within errors, there is little evidence for clumping on angular scales smaller than the VLA 
beam size of 50". Otherwise, one would expect a significant discrepancy between the T s calculated 
from the large VLA beam and the STIS pencil beam due to the filling factor of the gas. Finally, 
a visual examination of Figure 6 reveals no noticeable clumping in H I emission or H I velocity 
on angular scales slightly larger than 50". Thus, a reasonable estimate of the angular size of the 
absorbing clouds is ~ 50", which corresponds to a physical size of 5hj^ kpc at the distance of 
NGC 3067. This characteristic size yields a mean physical density for this cloud of ~ 0.007 cm -3 , 
remarkably similar to the nni = 0.004-0.006 cm~ 3 found from photoionization modeling by T99. 
This agreement is further evidence that the absorbing cloud is rather uniform on a 5h^Q kpc scale. 

It is surprising that a spherical cloud with a 5/iyQ 1 kpc radius would have a line width that is 
purely thermal with no contribution from turbulent or bulk motions. Galactic HVCs tend to have 
larger line widths than that found in H I 21 cm absorption for the 1420 km s _1 cloud; the 6-value 
for a typical Complex C sightline is ps 10 km s _1 for low ions like HI O I, Fe II, and S II (see 
Collins, Shull, &: Giroux 2003; Sembach et al. 2004). It is unclear why the H I 21 cm line width in 
the 1420 km s _1 cloud is so narrow (b = 2.5 ±0.1 km s _1 ). Perhaps the cloud is sheet-like rather 
than spherical such that it is smaller along the line of sight than in the transverse direction. 

Based upon new Ha emission detections, Putman et al. (2003) find that most Galactic HVCs 
lie within ~ 40 kpc of the Galaxy and are not extragalactic members of the Local Group. The 
Galactic HVC Complex C spans ~ 90° on the sky and lies at a distance of 2-14 kpc (depending 
on the spiral arm with which one assumes the HVC is associated; Putman et al. 2003), yielding a 
physical size of 3-22 kpc. This result implies that the absorbing clouds around NGC 3067 not only 
lie at comparable galactocentric distances as Galactic HVCs but also have comparable physical 
sizes. 
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3.3. Ionization Balance 

T99 found that all of the low and intermediate ionization species in the 3C 232 sightline could 
be explained by a photoionization model with an extragalactic ionizing flux of § lon > 2600 pho- 
tons cm~ 2 s _1 assuming a power-law spectrum with index 1.8 and no ionizing flux from NGC 3067 
(if one allows a significant contribution from NGC 3067 then it becomes impossible to constrain the 
extragalactic radiation field; T99). T99 also placed an upper limit on the escape fraction of ioniz- 
ing photons from NGC 3067 (f esc < 0.02) based on Ha images of the 3C 232/NGC 3067 system. 
This value is comparable to the escape fraction for Galactic HVCs (f eS c = 0.01-0.02) measured 
by Putman et al. (2003). The photoionization model of T99, which fits all of the low ions (Na I 
through Mg II and Fe II) predicts C IV and Si IV equivalent widths of ~ 10 mA, which would be 
undetectable in the G140M spectra presented here. Thus, the easily detected higher ions are much 
more likely to be collisionally-ionized by shocks or some other process despite sharing the same 
kinematics with the lower ions. 

In addition, the photoionization modeling of T99 sets a modest lower limit on metallicity for 
the ensemble of all three of these clouds of [Z/Zq] > —0.6 (Z > 0.25^0). This limit is specifically 
derived from the Si II 1304 A line (see Table 3) but is consistent with the observed strengths of 
other detected lines. This limit was derived assuming the extreme case of a purely extragalactic 
ionizing spectrum with an ionizing flux of \ v = 1 x 10~ 23 ergs s" 1 cm' 2 Hz -1 sr _1 as inferred by 
T99; however, if we assume the opposite extreme of an ionizing spectrum that comes purely from 
ionizing photons escaping the plane NGC 3067 with f esc = 0.02, the metallicity limit does not 
change appreciably. The metallicity limit for the NGC 3067 clouds of Z > 0.25^0 is comparable 
to recently measured metallicities of Galactic HVCs. Tripp et al. (2003) used a STIS echelle 
spectrum of the background quasar 3C 351 to find a metallicity of Z = 0.1-0.3 Z & for Complex 
C and Sembach et al. (2002) used a Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum of 
the background quasar Ton S120 to find a 3a metallicity limit of (O/H) < 0.46 solar for the 
compact high velocity cloud CHVC 224.0-83.4-197. Collins, Shull, k Giroux (2004) find similar 
metallicities of Z = 0.34 ± 0.14 Z & and Z = 0.06+g;^ Z Q for the highly ionized HVCs detected 
toward PKS 2155—304 at velocities of vlsr = —140 and —270 km s _1 , respectively. 

Highly ionized gas has been detected in several Galactic HVCs. Fox et al. (2004) have detected 
Complex C in a similar range of ionic species to those we detect for the clouds toward 3C 232. 
Using spectra from HST/STIS and FUSE, they found that O VI, Si II, and H I are all coincident 
in velocity, mimicking the ionization structure that we see in the NGC 3067 clouds. By modeling 
the high ion column density ratios iV(Si W)/N(0 VI), N(C W)/N(0 VI), N(N V)/N(0 VI), 
Fox et al. conclude that neither a single collisional ionization equilibrium nor photoionization 
equilibrium model can produce the highly ionized gas observed in Complex C. They hypothesize 
that the highly ionized gas is a collisionally ionized skin that forms in the conductive or turbulent 
interfaces that arise when the predominantly photoionized cloud moves through a surrounding hot 
medium. Given the similar ionization structure between Complex C and the NGC 3067 clouds, 
it is plausible that the same phenomena produce the ions in both; i.e., the low ions observed by 
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S91 and the intermediate ions observed by T99 have strengths and ratios consistent with a single 
photoionization model but not the Si IV and C IV, which require another ionization source. While 
Fox et al. discuss several alternative hypotheses for the observed ionization structure, the analogy 
between Complex C and these clouds is retained. We note that the inferred velocities of the three 
HVC-like clouds in NGC 3067 relative to a static halo in that galaxy yield shock temperatures 
easily sufficient to produce large amounts of C IV and Si IV. 

Lehner, Keenan, & Sembach (2001) have detected the Galactic HVC 291.2-41.2+80 along the 
line of sight to the Magellanic Bridge star DI 1388 in several species that we have detected toward 
3C 232, including C II, O I, Si II, Si III, Si IV, and S III. However HVC 291.2-41.2+80 has not been 
detected in H I to a 5a limit of 10 19 cm" 2 (Lehner et al. 2001), which suggests that this Galactic 
HVC is analogous to the 1370 and 1530 km s" 1 clouds around NGC 3067. Using the C Il/Si II 
and O I column densities to trace the hydrogen column of HVC 291.2—41.2+80 yields estimated 
column densities of log (iVm + A^hii) = 17.62 cm -2 and logiVHi = 16.24 cm -2 , respectively, which 
Lehner et al. argue indicates that much of the hydrogen is ionized since C II and Si II trace both 
neutral and singly- ionized gas whereas O I traces only neutral gas. These Galactic analogs suggest 
that the 1370 and 1530 km s~ x clouds could also be significantly ionized and not predominantly 
neutral as is the case for HVCs with higher H I column densities. 

The highly ionized HVCs detected by Sembach et al. (1999) toward the QSOs Mrk 509 and PKS 
2155-304 have a different ionization structure than the NGC 3067 clouds or HVC 291.2-41.2+80 
in that they are detected exclusively in high ions such as C IV, Si IV, and N V with no detectable 
absorption from low or intermediate ions (e.g., Si II or C II) or H I 21 cm emission at corresponding 
velocities. A search for 21 cm emission near these sightlines revealed H I HVCs at comparable 
velocities to the highly ionized clouds within ~ 2° of the sightlines (the H I HVC with the highest 
column density has 7V"hi < 2 x 10 18 cm~ 2 ; Sembach et al. 1999), which suggests that the highly 
ionized gas traces the extended ionized regions of the H I HVCs in much the same way as the highly 
ionized skin proposed by Fox et al. (2004) for Complex C. Sembach et al. also find that the lack of 
intermediate ions in their highly ionized HVCs argues against collisional ionization as the primary 
ionization mechanism for these clouds; instead they prefer a photoionization model in which the 
source possesses a hard QSO-like spectrum as the principle photoionization mechanism, although 
they don't rule out a contribution to the ionization balance from collisionally ionized gas. The 
requirement of a QSO-like source in their photoionization models implies that the highly ionized 
HVCs of Sembach et al. lie at large distances from the Galactic plane (30-200 kpc if the highly 
ionized HVCs and their H I counterparts are in thermal pressure equilibrium; Sembach et al. 1999) 
where they are exposed to the extragalactic radiation field. 

Subsequent observations of these HVCs with FUSE (Sembach et al. 2003) and STIS (Shull, 
Tumlinson, k, Giroux 2003; Kraemer et al. 2003) have changed the initial interpretation of the 
Sembach et al. (1999) GHRS data. Strong O VI absorption was found in the FUSE spectra of these 
HVCs at column densities that are too large for reasonable photionization models to accomodate. 
Many intermediate ions are detected in the STIS spectra, including Si II, Si III, C II, and C III; 
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H I was even detected in two HVCs toward PKS 2155-304 at N m = 10 15 to 10 16 cm" 2 (Collins 
et al. 2004). Collins et al. (2004) argue that the O VI in these HVCs is collisionally ionized and 
the presence of absorption in ions ranging from H I to C IV, Si IV, and O VI suggests that these 
HVCs are multiphase. They propose that the high ions in these clouds arise from hot gas that 
is potentially produced behind bow shocks from clumps falling toward the Galactic plane. This 
interpretation is remarkably similar to that of Fox et al. (2004) for Complex C and suggests a 
similar source for the high ions detected in the NGC 3067 clouds. 

4. Conclusions 

We have presented sensitive new radio and ultraviolet spectra to determine if the multiple 
velocity components (cz = 1317 ± 2, 1417 ± 2, 1530 ± 10 km s" 1 ) found in Na I and Ca II by S91 
and Mg I, Mg II, and Fe II by T99 toward 3C 232 are also present in H I Lya, Si IV, and C IV. Our 
H I 21 cm data from the VLBA and Arecibo show evidence for only the strongest (1420 km s _1 ) 
velocity component and set upper limits, Nm < 9 x 10 18 cm~ 2 , on the column density of the other 
two components. The 1420 km s -1 absorber is a Lyman limit system with a column density of 
A^hi = 10 20 ' cm~ 2 based on Voigt profile fitting of the Lya absorption feature in our STIS G140M 
spectra. Lyman limit systems are generally associated with the halos of luminous galaxies (Steidel 
1995, 1998; Stocke et al. 2005), and indeed our best value of the measured spin temperature for this 
absorber (T s = 435 ± 140 K) is significantly larger than the spin temperatures (~ 50 K; Akeson & 
Blitz 1999) found in the Galactic disk and comparable to the lower limit of T s > 180-320 K found 
toward one Galactic high velocity cloud (HVC; Akeson & Blitz 1999). The Na/Ca ratio for these 
absorbers is not consistent with gas in the Galactic disk either (S91), further suggesting that these 
absorbers are located in the halo, not the extended disk, of NGC 3067. 

While the 1420 km s _1 component is still the only component detected in H I, we did detect all 
three velocity components in both lines of the Si IV and C IV doublets in our STIS G140M spectra. 
The presence of all three velocity components in various metal species and ionization states implies 
that the highly ionized gas is spatially coincident with the neutral and singly ionized gas. We 
hypothesized that the three velocity components correspond to three clouds of gas ~ HhjQ kpc 
from the nucleus of NGC 3067 that are analogous to Galactic HVCs. We have found several lines 
of evidence that support this argument: 

• We find that galactocentric radial velocities of the NGC 3067 clouds (—95 to +65 km s" 1 ) lie 
in the velocity range expected for a Galactic HVC with vlsr = 90 km s _1 . 

• Assuming a cloud location along the minor axis of NGC 3067 and purely radial motion, the 
most likely kinematic solution for the 10 20 cm -2 cloud and the underlying disk of NGC 3067 
require this cloud to be falling toward the disk as frequently observed for Galactic HVCs. 

• Our estimated cloud size of 5/tyg 1 kpc is comparable to the estimated size of the Galactic HVC 
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Complex C. 

• The spin temperature of the NGC 3067 clouds (T s = 435 ± 140 K) is consistent with the lower 
limits found for Galactic HVCs (T s > 200 K; Akeson & Blitz 1999). 

• The photoionization modeling of T99 places a lower limit on the metallicity of the NGC 3067 
cloud ensemble of Z > 0.25 Z@, which is similar to the metallicities derived along several 
sightlines toward Galactic HVCs. 

• We detect the same intermediate ions in the NGC 3067 clouds (e.g., C II, Si II) that are 
detected in many highly ionized Galactic HVCs. 

• The upper limit on the escape fraction of UV radiation to these clouds is comparable to that 
found from the Galaxy to Galactic HVCs (Putman et al. 2003). A much larger escape fraction 
would be expected along the minor axis of a starburst wind (see Dove, Shull, & Ferrara 2000; 
MacLow & Ferrara 1999). 

Thus, the overall physical picture of the 3C 232/NGC 3067 system is that of a Galactic "fountain", 
not an outflowing starburst wind. 

We find no evidence that halo gas along this sightline is escaping the galaxy's gravitational 
potential, despite the fact that NGC 3067 is a modest starburst with a star formation rate of 
1.4 M yr" 1 . The 3C 232 sightline is located near the minor axis of NGC 3067, which is where 
one would expect to find signatures of an outflowing starburst wind (e.g., Heckman et al. 1990). 
Nevertheless, regardless of whether we assume that the halo gas shares the rotation direction of 
the disk, the deprojected velocities of the NGC 3067 clouds are all less than the escape velocity 
from the gravitational potential well if they are travelling on purely radial orbits along the minor 
axis of NGC 3067. Furthermore, the most consistent analysis of the disk/halo gas kinematics from 
a detailed H I emission map (CvG) determines that two of these clouds (including the H I 21 cm 
detected cloud) are infalling. For either kinematic solution (halo gas rotating or counter-rotating) 
the A^hi = 10 20 cm -2 system is gravitationally bound to NGC 3067. Thus, if these clouds truly are 
analogous to Galactic HVCs, one would expect that Galactic HVCs are also bound to the Galactic 
gravitational potential rather than that of the Local Group and that Galactic HVCs are located 
closer to the Galaxy than the Local Group barycenter. 

If instead we take the perspective that at least one of the velocity components observed in 
NGC 3067 towards 3C 232 is an outflowing starburst wind, do we observe what might be expected 
for such a circumstance? Almost certainly not. First, we might expect the higher ions to have 
different kinematics from the low ions; the latter would be expected to be seen at lower galactocen- 
tric velocities as cool clumps in a hotter, more highly ionized wind. While the kinematic solution 
does strictly allow an outflowing, even escaping wind if the geometry is just right (e.g., ejection 
nearly along the plane of the sky), the more straightforward kinematic solution is that two of the 
clouds are infalling, certainly not expected in a starburst wind scenario. The limit on the escape 
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fraction of ionizing radiation from NGC 3067 (/ esc < 2%) is also much lower than expected in a 
situation where an outflowing starburst wind would be expected to clear away almost all of the 
neutral material above the disk (Dove et al. 2000; MacLow & Ferrara 1999). The large inferred size 
of 5hjQ kpc for the neutral Nm = 10 20 cm~ 2 cloud is far larger than might be expected for a cool 
clump in a starburst wind; i.e., this size is nearly half of the cone opening angle (~ 65°) suggested 
by Heckman et al. (1990). Thus, a starburst model is a poor match to many of the observations 
for these clouds described in this paper. 

This result also suggests that the superwinds produced by other comparably massive local 
starbursts (e.g., NGC 253, NGC 4945) might not be capable of escaping their galaxy's potential 
well either. However, starburst winds may escape the shallower gravitational potential wells of less 
massive galaxies (L < 0.1 L*; Martin 2003; Stocke et al. 2004). Even though all of the low, inter- 
mediate, and high ions detected in these putative HVCs share the same kinematics, there remains 
the possibility that even more highly ionized gas could possess different kinematics that trace a 
wind escaping from NGC 3067. Unfortunately, at F A (1200 A) » 4-5 x 10" 15 ergs s" 1 cm" 2 A -1 , 
3C 232 is too faint in the far-UV to be observed easily with FUSE, which could be used to test this 
hypothesis by observing the O VI AA1032, 1038 absorption lines at high signal to noise ratio. 

The QSO/galaxy pair PKS 1327-206/ESO 1327-2041 shares many characteristics with the 
3C 232/NGC 3067 system. PKS 1327-206 and ESO 1327-2041 have a projected separation of 
13/ifo 1 kpc (CvG). ESO 1327-2041 has a similar luminosity to NGC 3067 (L » 1.2 L*; Lauberts & 
Valentijn 1989; Marzke et al. 1994) and was noted by Giraud (1986) and CvG to possess a polar ring 
of extraplanar gas. Thus, it is worthwhile to re-examine the gas kinematics and morphology of this 
QSO/galaxy pair in light of the HVC interpretation presented here for the similar 3C 232/NGC 3067 
system. 

The results we present here strongly support the model for Galactic HVCs that places them 
within a few kpc of the disk of the Milky Way (e.g., Bland-Hawthorn & Maloney 2002; Putman 
et al. 2003) and not at ~ 1 Mpc as "extragalactic clouds" within the Local Group (Blitz et al. 
1999). In the case of NGC 3067 there is little doubt that these absorbers are within a few kpc 
of its galactic plane and have remarkably similar properties to those of Galactic HVCs. Deep H I 
21 cm searches for HVC analogs in nearby galaxy groups (Verheijen et al. 2001; Zwaan 2001) and 
in non-targeted H I surveys (Zwaan k, Briggs 2000) have failed to find extragalactic H I clouds, 
which would be expected if the Blitz et al. (1999) hypothesis were correct. The present result adds 
to the significant evidence against the extragalactic evidence for HVCs. However, CvG note that 
the very extended H I in NGC 3067 and the possibly complex kinematics of the extra-planar H I 
are quite similar to the H I morphologies seen in interacting galaxy systems (e.g., M 51 &: Mrk 
348), although NGC 3067 has no obvious close companion. Could it be that the Milky Way and 
M 31 are also unusual in possessing such an extensive system of extra-planar H I clouds? 

The authors thank T. Ghosh and C. Salter for their help with the Arecibo observations and data 
reduction, T. Rector for his help with the VLBA reductions, and P. Maloney for helpful discussions. 
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